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Abstract 

Gravitational waves emitted by kinks on infinite strings are investigated using detailed estima- 
tions of the kink distribution on infinite strings. We find that gravitational waves from kinks can 
be detected by future pulsar timing experiments such as SKA for an appropriate value of the the 
string tension, if the typical size of string loops is much smaller than the horizon at their formation. 
Moreover, the gravitational wave spectrum depends on the thermal history of the Universe and 
hence it can be used as a probe into the early evolution of the Universe. 
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I. INTRODUCTION 



It is well known that cosmic strings are produced in the early Universe at the phase 
transition associated with spontaneous symmetry breaking in the Grand-Unified-Theories 
(GUTs) [l]. Although cosmic strings formed before inflation are diluted away, some GUTs 
predict series of phase transitions where symmetry such as U(1)b-l (B and L denotes 
baryon and lepton number) breaks at low energy and produces strings after inflation. Cosmic 
strings are also produced at the end of brane inflation in the framework of the superstring 



theory 



3j. Once produced, they survive until now and can leave observable signatures. 
Thus, the cosmic strings provide us with an opportunity to probe unified theories in particle 
physics which cannot be tested in terrestrial experiments. 

Various cosmological and astrophysical signals of cosmic strings have been intensively 
studied for decades. Especially, many authors have studied gravitational waves (GW) emit- 
ted from the cosmic string network, in particular, GWs from cosmic string loops. Cosmic 
stringloops oscillate by their tension and emit low frequency GWs corresponding to their 
size Q B. Moreover, string loops generica.lv have cusps and hinhs, and these structures 
cause high frequency modes of GWs, i.e. GW bursts [6H9(. 

On the other hand, GWs from infinite strings (long strings which lie across the Hubble 
horizon) have attracted much less interest than those from loops. This is because there 
are only a few infinite strings in one Hubble horizon, so GWs from them are much weaker 
than those from loops whose number density in the Hubble volume is much larger. 1 Loops 
can emit GWs with frequency u > (at)~ l , where a is the parameter which represents the 
typical loop size normalized by the horizon scale. Therefor, unless a is much less than 1, 
GWs from loops dominate in the wide range of frequency bands detected by GW detection 
experiments, e.g. pulsar timing or ground based or space-borne GW detectors. 



However, there are several reasons to consider GWs o 
papers suggest that a is much less than previously thought 



infinite strings. First, recent 
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llj , and hence loops cannot 



emit low frequency GWs which can be detected by pulsar timing experiments. If this is 



1 For cosmic strings whose reconnection probability is small, such as cosmic superstrings, many infinite 
strings can exist in a Hubble horizon. For such kinds of cosmic strings, however, the number of loops 
in a horizon is accordingly large, hence the situation that the loop contribution is stronger than infinite 
strings remains unchanged. 
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true, GWs from infinite strings can give dominant contribution to the GW background at 
those low frequencies, as we will show in this paper. Second, infinite strings can emit GWs 
with wavelength of horizon scale because the characteristic scale of infinite strings is as long 
as the size of the horizon. This may be detected by future and on-going CMB surveys. 
Loops cannot emit GWs with such long wavelength. 2 Thus, we are lead to consider GWs 
from infinite strings. 

A structure on infinite strings which is responsible for GW emission is a kink. 3 Kinks are 
produced when infinite strings reconnect and kinks on the infinite strings can produce GW 
bursts. The kinks are sharp when they first appear, but are gradually smoothened as the 
Universe expands. Moreover, when an infinite string self-intercommutes and produces a loop, 
some kinks immigrate from the infinite string to the loop, and hence the number of kinks 
on the infinite string decreases. Recently, Copeland and Kibble derived the distribution 



function of kinks on an infinite string taking into account these effects 13( . However, they 
did not study the GWs from kinks. 

Therefore, in this paper, we investigate GWs from kinks using the kink distribution 



obtained in |13[ and discuss the detectability of them by future GW experiments, especially 
by pulsar timing, such as Square Kilometer Array (SKA) and space-based detectors such as 
DECIGO and BBO. 

This paper is organized as follows. In section 2, we briefly review a part of the basis 
of cosmic strings. In section 3, we derive the distribution function of kinks on infinite 
strings according to In section 4, we derive the formula of the energy radiated from 

kinks per unit time. In section 5, we calculate the spectrum of stochastic GW background 
originating from kinks, using formula derived in section 4. Section 6 is devoted to summary 
and discussion. 



2 Other recent simulations imply a ~ 0.1, which is much bigger than G/j,. In this case loops can emit GWs 
whose wavelength is comparable to or somewhat shorter than the horizon scale [ljj ■ 

3 Generically, infinite strings do not have any cusp since its existence depends on the boundary condition 
of the string, i.e., the periodic condition. 
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II. DYNAMICS OF COSMIC STRINGS 



The dynamics of a cosmic string, whose width can be neglected, is described by the 
Nambu-Goto action, 

S = -ft J rf 2 Cv / -det( 7afe ). (1) 

where ( a (a = 0, 1) are coordinates on the world sheet of the cosmic string, 7^ = g iiv x >i a x v h 
( xf a = ) i s the induced metric on the world sheet, and fi is the tension of the string. The 
energy-momentum tensor is 

T"( X ) =nj d 2 C^-det{ lab)l ab x^ b 5\x - X(()), (2) 

where X = X(() is embedding of the world sheet on the background metric. If the back- 
ground space-time is Minkowski one, we can select the coordinate system (J 4 ) = (r, a) 
which satisfies the gauge conditions 

r = t (physical time), x r • x CT = 0, x\ + x\ = 0. (3) 

The time scale of a GW burst is much shorter than the Hubble expansion, hence we consider 
an individual burst event on the Minkowskian background. The general solution of the 
equation of motion derived from the action (CQ) is 

x? = ^(a"(u) + a' 2 (u) = b'» = 1 (4) 

where u = a + t,v = a — t. We call a(u) (b(v)) the left (right)-moving mode. Then, Eq. (j2J) 
can be rewritten in terms of a(u) and b(v), 

T^(k) = ^m)r_{k) + imim), (5) 

I1(k) = j dua'^(u)e ik - a{u)/2 , It{k) = J dvb'^vy*^ 2 , (6) 
where T» v {k) is the Fourier transform of the T^(x), i.e. T» v {k) = J d 4 xT^(x)e ihx . 

III. DISTRIBUTION FUNCTION OF KINKS 

Kinks can be defined as discontinuities of a' or b'. They are produced when two infinite 
strings collide and reconnect because a' and b' on the new infinite string are created by 
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connecting a's or b's on two different strings. Let us suppose that a' jumps from a'_ to a' + 
at a kink. Then the "sharpness" of a kink is defined by 



^ = -(l-a' + -a'_). 



(7) 



Thus the norm of the difference between a'_ and a' + is | Aa'| = 1\fip. The production rate 
of kinks is given by 13] 



production 



AV 



gW, 



(8) 



where N(t, ip)dip denotes the number of kinks with sharpness between ip and ip + dip in the 
volume V, A and 7 are constants related to string networks, whose values are [yj], 



A r ~ 0.20, A m ~ 0.21, 7r ~ 0.3, 7 m ~ 0.55, 
Here the subscript r(m) denotes the value in the radiation(matter)-dominated era. 

g(iP) = —V^(l5-6^-iP 2 ) 

25o 



(9) 



(10) 



and we set g(ip) = for ip < 0, 1 < ip. The correlation length of the cosmic strings ^ is given 
by £ - it. 

Produced kinks are blunted by the expansion of the Universe. The blunting rate of the 
kink with the sharpness ip is given by 13] 



t 
iP 



-2(r\ 



(11) 



stretch 



where ( is a constant which, in the radiation(matter)-dominated era, is given by ( r ~ 
0.09 (C m ^ 0.2). 

On the other hand, the number of kinks on an infinite string decreases when it self- 
intercommutes since some kinks are taken away by the loop produced. The decrease rate of 



kinks due to this effect is given by 13] 



N 
N 



V 



(12) 



to loop 

It' 

where 77 is constant which, in the radiation(matter)-dominated era, is given by rj r ~ 
0.18 (fA»~0.1). 

Taking into account these effects, the evolution of the kink number iV obeys the following 
equation, 

N = ^rrgiip) + —-^-(ipN) — —N. (13) 





definition 


radiation dom. 


matter dom. 
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0.3 


0.55 


c 








0.09 


0.2 


d 








1.1 


1.2 


A 




2C-/3 




-0.92 


-0.8 


B 




4C-/3 




-0.74 


-0.4 


C 




8C)/[3(jS - 


20] 


0.14 


-0.17 



TABLE I: Various constants appearing in the calculation of kink distribution and the spectrum 
of GWs from kinks. 



We have to solve this equation under an appropriate initial condition. The initial con- 
dition to be imposed depends on how strings emerge. We consider two typical scenarios. 
As a first scenario, let us suppose that cosmic strings form when spontaneous symmetry 
breaking (SSB) occurs in the radiation dominated Universe. After the formation, cosmic 
strings interact with particles in thermal bath, which acts as friction on the string. At first 
this friction effect overcomes the Hubble expansion. (This epoch is called friction-dominated 
era.) The important observation is that the friction effect smoothens strings and washes out 

n 

the small scale structure on them [14{. The temperature at which friction domination ends 



is given by 



15| 



T c ~ G>M pl , (14) 



where G is the Newton constant and M p i is the Planck scale. If the temperature at the 
string formation is lower than this critical value, the friction effect can be neglected from 
the formation epoch to the present day, and kinks emerge at the same time as appearance 
of strings. 

As a second scenario, let us suppose that cosmic strings are generated at the end of 
inflation, like strings made by condensation of waterfall fields in supersymmetric hybrid 



inflation 



16 



18], and cosmic superstrings left after annihilation of the D-brane and anti 



D-brane in brane inflation 2|, |3[. In this case, cosmic strings form in the inflaton-oscillation 
dominated era, which resembles the matter-dominated era. Although strings feel frictions 
from dilute plasma existing before the reheating completes, the effect is negligible as long 
as the reheating temperature after inflation T r is lower than ~ GfiM p \. In this case, kinks 
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begin to be formed right after the formation of strings. Otherwise, the temperature at which 
kinks begin to appear is given by Eq. (I14p . 

In anyway, even if such kinks survive friction, they become extremely dense in the later 
period, so that they do not affect the observable part of the GW spectrum, as we will see 
later. On the other hand, if the reheating temperature is lower than T c , the kinks from the 
first matter era definitely continue to exist, and if the reheating temperature is extremely 
low, GWs from these kinks can be observed. Therefore, we concentrate on the situation 
where the reheating temperature is low enough so that the kinks generated in the first 
matter era survive without experiencing the friction domination. We denote the time when 
kinks starts to be formed by 

Then we can get the solution, but its precise form is very complicated. Here we assume 
that the matter-dominated epoch follows after inflation, and the reheating completes at 
t = t r after which the radiation dominated era begins. If we focus on only the dominant 
term and neglect 0(1) numerical factor, we get 




for ip > (^) 
for ip < (tf) 



(15) 



in the first matter era, 




for ^>><0i (t) 

for ^f D) (t) <^<^f D) (t) 



(16) 




in the radiation era where 




(18) 



(17) 
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and 



dN , M 



^-p m /2U t -i f or ^,>^ MD )(t) 

/^m £r / Cm 



L)^-^ 4 - ^-W r i for ^» (f) < ^ < 

r M2U t -i for D )(t) < ^ < V»f D) (t) 

, / \ Pm~\~Qm / i \ Pr~\~C,r / , \ Pm~\~Cm f\JT\\ 

^ (x) (£) (t) r* for ^ < ^f D) (t) 

(19) 



in the second matter era, where 



2Cr, 



^r\t) = (^j , (20) 

/ + \ 2 Cm / . \ 2C r 

# D) w = m fM , (21) 



t J \ t eg 



t\ 2^777, / , \ 2£ r / , ^ 2c^ 
PCI \ I "V \ I " 



Here f3 is the constant related to the string network (/3 r ~ 1.1, /3 m ~ 1.2), and t eq denotes 
the matter-radiation equality epoch. We have converted N, distribution in the volume V, 
to dN/dip = N(t,ip)/(V/^ 2 ), which is the distribution per unit length. The derivation of 
the above expression of dN/dip is described in Appendix LAI 

Let us consider the physical meaning of this distribution function. It is not difficult to 
consider how the number of kinks in the horizon changes as time goes on. When a kink is 
born, its sharpness ranges from to 1, but the typical value is 0(0.1). Therefore, the kink 
distribution has a peak at ip ~ 0(0.1) at the very early stage. Then kinks are made blunt 
by the cosmic expansion. Thus newly produced kinks are sharp (large ip) and old ones are 
blunt (small ip). Therefore, the peak value of ip is much smaller than 0.1 at the late stage. 
This peak consists of the oldest kinks. Fig. [T] roughly sketches the shape of the distribution 
function in the second matter era. 

It should be noted that the above distribution is derived without considering gravitational 
backreaction. Since the most abundant kinks are extremely blunt, they might be influenced 
by the gravitational backreaction and disappear. However, it is difficult and beyond the 
scope of this paper to take backreaction into account. We will comment on this issue in 
Appendix O 
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FIG. 1: The distribution function of kinks on infinite strings produced at the end of the inflation 
in the second matter era. i/j^ D \ , are given by Eqs. ([20]) - ([22]) . 

IV. THE SPECTRUM OF GRAVITATIONAL WAVES FROM KINKS 



Given the energy-momentum tensor of the source, one can calculate the energy of the 
GW in the direction of k with frequency u as [ljj] 



(k) = 2G / dcuu 2 T^{k)T^{k) - -\mk)\ 



dE 

POO 

= 2GK ijXm $) / duu 2 T ij *(k)T lm (k) (23) 
Jo 

1 „ „ . . l l . . 1 

A-ij,im{k) = SuSjm — 2kjk m 5u + —kikjk[k m — —Sij5i m + —Sijkik m + —Si m kikj (24) 

Thus by substituting the energy- momentum tensor dSJ), we find the energy of the GW 
after computing integral I± in Eq. (jSJ). It is given by 



dudVt 



F* 1 2 



\i+\ 2 \k-L\ 2 -\L\ 2 \k-f + \ 2 



\i+-r , 

k ■ I+\ 2 \k ■ I-\ 2 

- (i + ■ r*_)(k ■ r + )(k ■ j_) - (i; • /_)(£ • i + )(k • r 
+ (/+ ■ /_)(£ ■ • ii) + (/; ■ r )(k ■ f + )(k ■ L 



(25) 
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where /± = (4(*),4(*)>4(*))- 

The method to calculate I± is described in [7|, |20|. In the limit u -> oo, I± exponentially 

reduces to 0, unless at least one of the following conditions on the integrand is met. One is 

the existence of discontinuities of a' l (u) (or b n (v)). The contribution of a discontinuity of a' 

to I+{k) is 



a'i a H 



J* (Jfe) ~ ± ^ e M«.-*-a.j/2 (26) 

+ iw \1- fc-a^ l-Jfe-ai ' 

where is the position of the discontinuity, a* = a(it*) and we assume a' jumps from a' to 

a' + at m = u*. The region of length ~ u^ 1 around u = u* contributes to this value. We find 

~ ip/w, where ip denotes sharpness of the kink. The other condition is the existence of 

stationary points of the phase of the integrand, i.e. uj(u — k-a(u))/2 (or u(—v — k-b(v))/2). 

This condition is expressed as 

1 - jfe • a'«) = (or - 1 - k ■ h'(v s ) = 0) (27) 

at the point u = u s (v = v s ). The contribution of the stationary point of the phase to I+(k) 
is 

P + (k) ~ 4r 3 <(us)e^ u °-~ k - aM)/2 . I -4r(2/3). (28) 

w 1 \\k ■ a'" (ii s ) | / v <j 

For any value of a' (b'), there is one direction k that satisfies (1271) . i.e. k — a! (— b'). 
Therefore, every point of u can contribute to P + (k) for one direction k. 

One can find the energy of the GW burst from ONE kink by picking up the contribution 
from the discontinuity for one of I± (say, I + ) and from the stationary point for the other 
(say, 7_). The energy emitted when the kink is located at the world sheet coordinate 
(u,v) = (u*,v a ) is evaluated by substituting ( 126]) into I + and ( 1281) into J_ in ( |25i) . Then we 
obtain 



dE < KM r2 (l)^ 2 -i / 12 

(w, -wb'J = — ^-GiTu a ' 



dlMw v 3 Vlb'-b 



2a' • a' 1 



+ 7^—^ b 



//a 



[l + K-s! + Y (l + b' s -a' + )(l + b' s -a'_) (l + b' s -a'_)^ 

(29) 

where the subscript s represents the value at v = v s . The energy is radiated at every 
moment toward the direction of — b' from the kink. This is the formula which represents 
energy emitted in a short period in a small solid angle. The total energy emitted in the 
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short period At is found by multiplying An, which is the solid angle that the GW sweeps in 
this short period. The extent of the radiation has the solid angle ~ ((cu/lb"!) -1 / 3 ) 2 Q. The 
variation of the direction of the GW is roughly estimated by |Ab'| ~ |b"At|. As a result, 
we obtain 

An ~ (w/IKD-^lb^A* (30) 
and the energy emitted per unit time is 



duj 3 \ lb' c • b 

x 



2a' -a' 1 \ (31) 

Za ^ a - _i I I 1 10/3 



1 + h' s ■ a' + ) 2 (1 + b' s ■ a' + )(l + b' s • a'_) (1 + b' s ■ a'_) 2 

The terms in the large parenthesis in the 2nd line of Eq. (I3~lj) can be estimated by tak- 
ing average over the angle between the left- and right-moving mode as ( ( , 1+b , 1 a , ~ 
( (i+b' -a' )(i+h' a' ) ) ~ ^- magnitudes of b" and b' s " should be t~ l ) and £" 2 (~ t~ 2 ). A 

After making these substitutions and neglecting 0(1) numerical factors, we find 

lOG^cu-^t- 2 / 3 . (32) 



dP 

duj 



one kink 



So far we have evaluated the GW spectrum from one kink. However, there exist many 
kinks on an infinite string and the final observable GW spectrum is made from sum of 
contribution from these kinks. Therefore, /+ picks contributions of many kinks. Formally, 

4(*)=E4,m(*), (33) 
m 

where an integer m labels each kink. Thus Eq. (J25]) has cross terms of the contributions from 
different kinks, e.g., 7 +jm • I+ n - However, such cross terms must vanish since a GW burst 
from a kink is a local phenomenon which relates only the region around the kink. In fact, 
the structure of the string around a specific point arises as a result of nonlinear evolution of 
the string network, and hence the values of b", b' s " and so on, are stochastic. We show that 
ensemble averages of cross terms vanish as expected, i.e. (P + m I+ n ) = in Appendix IB] 
where it is also shown that the kinks that dominantly contribute to the power of GWs with 



4 Although b' has small and dense discontinuities which represent kinks, v s is generally not a kink point. 
So the rough shape of h' s is determined by the global appearance of the string network and hence the 
length scale of its variation is roughly the curvature radius of the network. 
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frequencies ~ u are ones which satisfy 



dNy 1 



(34) 



In other words, if the interval of kinks with sharpness ~ ip is similar to the period of the 
GW under consideration, these kinks make dominant contribution to the GW. In the first 
matter era, using Eq. ( TT5"j) , Eq. (I3"4"|) simplifies to 5 

i> ~ (oot) 2 ^ /Am (35) 

for w < (tJt) Am t~ l . (Here, we set A = 2£ - A r ~ -0.92, A m ~ -0.8.) For w > 
(t*/t) Am t _1 , fV'f^) > * s satisfied by an arbitrary value of if). In Appendix IB1 it 
is shown that for u > (t^/t)" 4 " 1 ^" 1 the main contribution to comes from the kinks 
corresponding to the peak of the distribution, i.e. if) ~ {t*/t) m - If we denote the value of 
sharpness of kinks which make dominant contribution to as 4> max (tu,t), it is given by 

( Ut yc m /A m for u < ^i^A mt -i 

(^) 2Cm for U) > (u/ty^t- 1 

in the first matter era. In the radiation era, if) ma , x (u),t) is found in a similar way as 

\ ut yCr/A r for W<a; ( RD )( t ) 



^max(w,t) 



(36) 



(37) 



Ut) 



2C„ 



for u>ui RD \t) 



where 



to 



(RD) 



(RD) 



(t) 



(t) 



In the second matter era, ^ max is estimated as 



^max(w,f) ~ < 



^ t yCm/A m 



2D/A„ 



for W <^ MD) (t) 

for ^ MD) (t)<u;<4 MD) (t) 

for a;f D) (t)<^<4 MD) (t) 



(38) 
(39) 



(40) 



for > uf 1 D 



Strictly speaking, Eq. has two solution for ip, but it is sufficient to take larger one. 
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where 



(MD) 



U < MD) (t)=(^j [± 



uj 



t 



eq 



(41) 
(42) 
(43) 



Here, we set D = (3 r ( m — (3 m ( r ~ 0.11. 

As a result, assuming that the powers of GW from different kinks are roughly same as 
far as their sharpnesses are in the same order, (in other words, assuming that the quantities 
concerned with kinks, such as b", except their sharpness, are roughly same) we can estimate 
the total power of GWs with frequencies ~ u from all of the kinks in a horizon as, 



dP 



~ 00- 



tot 



dP 

doj 



i , , \ \ ,dN 

(W,^max(w,t)) X 1p — 
one kink t 



X t. 



(44) 



(45) 



The first factor denotes the power of GWs from one kink. The second factor denotes the 
number of kinks which satisfy ip ~ ^ max (w, t) per unit length. The third factor is length of 
an infinite string in a horizon. Then we finally obtain 

dP J 10G/i 2 M) Cm for r 1 < w < {U/t) A ™t~ l 

duJ tot |iOG/i 2 (if) 3 " 1 (tut)- 2 / 3 for u > [U/t^t- 1 
in the first matter era, 

' 10Gfi 2 (iot) c r for r 1 < w < o;f D) (t) 

lOG/i 2 (k) 2 ^/^ ( wf )CU for w CBD) (t) < ^ < w (RD) (t) 

^10G/i 2 (|) i?m (^) Br (a;t)- 2 / 3 for w > D) (t) 
in the radiation era, and 



dP 

■ 

doj 



tot 



(46) 



a; 



da; 



10G/x 2 (a;t) c '™ 



tot 



\ -2D/A m 

10G/i 2 (^) (a,t) c " 



Br 



for r 1 <a;<a;! MD) (t) 

for a,j MD) (t) < a; < a;f D) (t) 

for a;f D) (t)<a;<a;f D) (t) 



(47) 





v*0 




) 





a; 



(MD) , 



'3 'W 

in the second matter era, where B = A( — (3 (B r ~ —0.74, B m ~ —0.4) and C = (/3 — 
80/3(0 - 2C) (Cr - 0.14, C m ~ -0.17). A robust lower bound on the frequency of GWs is 
set to be t _1 corresponding to the horizon scale. 
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V. THE STOCHASTIC BACKGROUND OF GRAVITATIONAL WAVES FROM 
KINKS 



We have found Eqs. (j4"5]) . (l4l)]) and (117)1 as the total energy radiated per unit time in a 
horizon from kinks on an infinite string. Now we are in a position to calculate the density 
parameter of GWs defined by 

where p c denotes the critical energy density of the present Universe. Noting that the energy 
density of GWs decreases as a~ A and the frequency redshifts as a -1 , we get 



u'=u)Xao/a(t) 



a(t) 
a 



(49) 



where a(t) is the scale factor and ao represents its present value. Using Eqs. ( I4"5]) .( l4"6]) and 
( )47|) . we get the spectrum as 



f2g w (k->) 



607r(G/i) 2 (wt ) Cm for t^ 1 < ou < u ± 

-3D/A r 

-4D/A m / >. 4D/A m / \ -4D/A, 



60n(Gfi) 2 ( gt. ) (cot ) Cr for wi < u < cu 2 



6M<W (fe)" " (^) " (^)~ " (a*,)*- for u, 2 < u, < c 3 

, / ^ \ 3B m /2— 4_B r / \ 2B r / \ 2B m —2B r / \ — _B m 

^(G^W^ (fe) (4) (tO (fe) fOT ->^ 

(50) 



fir 



where 



! MD) (io)=(^) (J 1 , (51) 



t) (£) (£) i »' (52) 



(O \ 3A m /2— 2A,. / rp \ 2/t r / rji \ 2A m —2A r / tt \ —Am 
a (4) ® © * (m) 

where to and To denote the present age and temperature of the Universe, T r is the reheating 
temperature and is the Hubble parameter at the end of inflation. 6 This formula has 



6 Here we neglected the effect of cosmological constant. The GW spectrum will be slightly modified if 
the cosmological constant is taken into account. Detailed estimation of this effect is beyond the scope 
of this paper since it needs a simulation of cosmic string network evolution in the cosmological constant 
dominated Universe. 
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FIG. 2: A schematic picture of gravitational wave spectrum from kinks (|50p . uj\, 0J2, ^3 are given 
by Eqs. (j5B - ([551 h 



complicated exponents. Using (3 and ( with values given above, Q m /Q r ~ 5.5 x 10 3 and 
T /M pl ~ 9.6 x 1(T 32 , Eq. §0§ simplifies to 



10 2 (G/i) 2 (a;t )- a17 
10(G» 2 (u;i ) - 14 



where 



10 18 (G/i) 



. 0.56 

2 I T r \ I, . \-0.17 



for 1 < u < 0J\ 
for u\ < u < UJ2 
for U2 < co < W3 



lO^G/i) 2 !^) ' 68 ^) ' 4 ^)" 2 / 3 for a;>^ 3 , 



(54) 



0J\ ~ 

CJ3 ~ 



10V. 

ID- I I 



l.S 



t: 



10 



51 



T,. 



M 



pi 



-1 

> 

0.24 



M, 



pi 



0.8 



(55) 
(56) 

(57) 



We find that the integral in Eq. ( 149]) is dominated by the contribution from the period 
near the present for the whole range of u. In other words, almost all of the present energy 
of GWs from the kinks on infinite strings comes from those radiated around the present 
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epoch. This does not mean that the kinks produced at present make dominant contribution 
to the energy of GWs. We see above that for given frequency u, the kinks which dominantly 



21|, the second inflection point falls in the 



contribute to u£ are determined by Eqs. (I3"6"|) . (l3"7j) or (HUj) . Accordingly high frequency 
modes arise from dense, blunt and old kinks, and low frequency modes arise from thin, 
sharp and new ones. The first line of Eq. (150j) corresponds to GWs from new kinks which 
were born after the matter-radiation equality, the second line corresponds to GWs from 
old kinks which were born in the radiation era, the third line corresponds to GWs from 
older kinks produced between the end of the inflation and the start of the radiation era and 
the last line corresponds to GWs which came from the most abundant and oldest kinks. 
Figure [2] sketches the shape of the spectrum. The spectrum has three inflection points at 
U\,U2 and U3. This is due to change of the type of kinks which mainly contribute to f2 gw . 
Positions of these inflection points depend on the reheating temperature T r . If T r takes 
the value around its lower bound, say, 10 MeV 
observable region, as we will see. Even if T r is so low, the natural value of ii* makes the 
third inflection far above the observable region. If the reheating ends immediately and ii* 
is as small as possible, the region between the second and third inflection is so short that 
the third inflection enters the observable region. 

This is a crude estimation, and in order to derive the realistic spectrum we should take 
into account subtlety described in [?], ls| where the authors claimed that GWs from kinks 
are burst-like and hence GW bursts with rare event rate ("isolated" GWs ) should not be 
counted as constituent of the stochastic GW background. We should calculate the stochastic 
GW background spectrum using the following formulae [z, 8] 

*W/) ~ —(ft ) 2 h 2 coni (f) (58) 

h 2 con( (f) = J tj9W, z) - l)n(/, z)h 2 (f, z) (59) 

= 7^ (60) 

dN ~ ^9^,z)(l + zy^ max (u z ,t)N(ij max (u z ,t),z)t-\z)dV(z), (61) 
6 m (f,z) = [(l + z)ft(z)]-V s , (63) 
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dV 



547rtjjw§^(l + z)- 9 l 2 dz (1st matter era) 

72ntl^J^(l + z)~ 5 dz (radiation era) (64) 

547r*g((l + z) 1/2 - 1) 2 (1 + z)~ ll / 2 dz (2nd matter era) 



t(z) 



W^O- + 2 )~ 3/2t o (1st matter era) 
[1 + z)~ 2 t (radiation era) 

(65) 

'1 + z)~ 3 / 2 t (2nd matter era) 

where f = u/2ir and u z = u(l + z). dV means the proper spatial volume between redshifts 
z and z + dz. t(z) represents the cosmic time at the redshift z. n(u, z) represents the 
number of GW bursts at redshift ~ z with frequencies u superposed in a period of ~ 
N(ip, z)dVdip is the number of kinks with sharpness ip ~ ip + dip in the volume dV at redshift 
z, so N(ip, z)t 2 (z) ~ ^(t(z),ip). Isolated GW bursts are excluded from the calculation by 
inserting the step function in the integral in Eq. (159]) . h(f,z) is the logarithmic Fourier 
component of the waveform of the GW burst from one kink located at redshift z and can 
contribute to GW with frequency /. 

The results of calculation are shown in Figure[3]and FigurelU We take Gfi ~ 10~ 7 , close to 



the current upper bound from CMB observation [22j, and show the spectrum with frequency 
uj from the band of CMB experiments to that of ground-based GW detectors. These figures 
include both our crude estimate fl50l) and the improved one given by Eqs. fl58|) - fl65l) . We see 
that the latter is much smaller than the former. This is because the spectrum is dominated 
by GWs emitted recently, and recent GW bursts have a more tendency to be isolated. The 
fact that the difference between the two estimates becomes larger in higher frequency band 
might disagree with intuition, since the kinks corresponding to high frequency GWs are 
more abundant. However, the higher frequencies of GW are, the smaller the possibility that 
they overlap, because the period of oscillation becomes shorter and the extent of the GW 
beam becomes narrower. As a result, higher frequency GWs are more likely isolated in time. 

In Figure [3] we assume that cosmic strings were born by SSB at GUT scale and first 
kinks appeared at the temperature 10 12 GeV (JHJ). In Figure |4] we assume that strings 
were produced at the end of the inflation with extremely low reheating temperature (T r = 
10 MeV). This assumption makes the second inflection point of the spectrum visible in the 
observable frequency band. We also assume the inflation energy scale is sufficiently high so 
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FIG. 3: f2 gw in the case where strings emerge at the phase transition in the radiation era for 
G/U = 10~ 7 and ~ 10 12 GeV. The upper line represents the estimate using Eq. (|50|) . i.e. including 
"rare bursts", and the lower line represents the estimate using Eqs. ([58]) - ([63]) . i.e. excluding "rare 



bursts" . Sensitivity curves of various experiments are shown. That of D 



That of BBO correlated is derived from 



261 ] . Others are derived from 



1CIGO is derived from 



that third inflection point on the spectrum is far from the observable region. If we could 
observe the second inflection, we can deduce the reheating temperature. 

Note that the spectrum depends on G[i via the overall factor (Gfi) 2 . Therefore, when we 
vary the value on G/j,, the spectrum only moves upward or downward, and its shape (e.g., 
the position of bending) does not change. This situation is different from the case of GWs 
from cosmic string loops. In the case of cosmic string loops, since different values of G/i give 
different values of lifetime of loops, the resultant spectral shape is different % la]- 

Let us discuss detectability of this GW background. We have to see whether GWs from 
kinks exceeds not only thresholds of various experiments but also GWs from other source. 



The spectrum of GWs from loops was discussed in 



4ri9|, and the contribution of loops to 
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frequency[Hz] 



FIG. 4: figw in the case where strings emerge at the end of inflation for G\i = 10~ 7 , T r ~ 10 MeV. 
The upper line represents the estimate including "rare bursts", and the lower line represents the 
estimate excluding "rare bursts". 



GW background is much larger than that of kinks if they coexist in some frequency band. 
However, a loop cannot emit GWs with frequencies smaller than the inverse of its size. Thus 
there is a cut on the low frequency side of the spectrum of GWs from loops corresponding 
to the inverse of the loop size ~ (at) -1 . If a ~ Gfi, the spectrum of GWs from loops begins 
to appear at u > 10~ 12 Hz, and this covers the frequency band where both pulsar timing 
arrays and GW detectors have good sensitivity. However, a is one of the most unknown 



12], a may 



parameters in the cosmic string model. According to some recent simulations 
be much greater and the broader region may be covered by loops' GW. On the other hand, 
some recent studies 



10 



11] show the possibility that a is extremely small, say, a ~ {Gfi) n 



with n > 1. In such a case, GWs from loops dominate only very high-frequency region and 
GWs from kinks may be observable at low-frequency region. For example, if a < 10~ 9 , the 
band of SKA 23] can be used for detection of GWs from kinks, and for Gfi 



10 . g-^y 
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reaches the sensitivity of SKA. In a more extreme case a < 10 , the sensitivity band of 



space-borne detectors, BBO 24J and DECIGO 25j are open for detection of kink-induced 
GWs, and f2 gw exceeds the sensitivity of correlated analysis of BBO for Gfi ~ 10~ 7 , and 
that of ultimate-DECIGO for G[i > 10 _9 . 7 In even more extreme case, a < 10 -16 , it may 
be possible to detect the inflection point in Figure 4 for Gfi > 10~ 10 and determine the 
reheating temperature as T r ~ 10 MeV. 

Moreover, GWs from kinks may be detected through CMB observations. As opposed to 
GWs from loops, kinks can emit GWs with wavelength comparable to the horizon scale. 
These GWs induce B-mode polarizations, which is a target of on-going and future CMB 
surveys. The spectrum of GWs from kinks is quite different from inflationary GWs and 
hence its effect on CMB is also expected to be distinguished from that of inflationary origin. 
We will study this issue elsewhere. 

VI. CONCLUSIONS 

In this paper, we have considered gravitational waves emitted by kinks on infinite cosmic 
strings. We have calculated the spectrum of the stochastic background of such gravitational 
waves and discussed their detectability by pulsar timing experiments and space-borne detec- 
tors. It is found that if the size of cosmic string loops is much smaller than that of Hubble 
horizon, some frequency bands are open for detection of GWs originating from kinks. It 
can be detected by pulsar timing experiments for Gfi > 10~ 7 , and by space-borne gravita- 
tional wave detectors for much smaller Gfi, although the latter may be hidden by the loop 
contribution unless the typical loop size is extremely small. If it is detected, it will pro- 
vide information on the physics of the early Universe, such as phase transition and inflation 
models. Moreover, the spectrum shape depends on the thermal history of the Universe, and 
hence GWs from cosmic strings can be used as a direct probe into the early evolution of the 
Universe. Notice that the inflationary GWs also carry information on the thermal history 
of the Universe 
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31] . Although the inflationary GWs are completely hindered by GWs 
from cosmic strings if the value of G/i is sizable and GWs from kinks come to dominate in 
low-frequency region, GWs from cosmic strings also have rich information on the physics of 



7 Note that stochastic GWs from astrophysical sources such as white-dwarf binaries make a dominant 



contribution for to < 1 Hz 2Sj and the observable frequency range is somewhat limited. 
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the early Universe. 



Appendix A 

Here we derive the expression of dN/difi given in section 3. First, we change the variable 
from -0 to tpi = ijj(t/ti) 2< * in Eq. ( TLB"]) , where ti is the time when we set the initial condition. 
Then Eq. ( TIB"]) becomes 

tNtyi, t) + - 2Cj Ntyi, t) = ( (^j C A , (Al) 

where the dot now denotes the time derivative at constant ipi. This equation can be easily 
integrated to obtain 

N(M _ A f dt> \ (tA^ Ndt/t^u) (A9] 

For the distribution function during the first matter era, we set the initial condition at U = t* 
as N(ip,Q = 0. Then Eq. (TA2]) becomes 

N±M Am [* dt> (ft\ 2<m ,\ 

v(t) ^W/«-o^Hw ' 



By substituting Eq. ([TO]) into g in Eq. (IA3j) . performing the integration and omitting terms 
except for dominant one, we get Eq. ( TT5l) . For N/V during the radiation era, we set U = t r 
and get 

t) _ A r f dt< _ g ( / A ^ \ + /t,\ N ((t/t r )^,tr) ^ (A4) 



We use Eq. ( fl5l) for N(ip,t r ). Then Eq. ( 1A4I) simplifies to Eq. ( fT6|) by picking only the 
dominant term. The expression during the second matter era (Tl9l) can be obtained in the 
same way. 

Appendix B 

Here, we prove kinks which dominantly contribute to GWs with frequency u are those 
which satisfy Eq. (|34|) . and evaluate the integral in Eq. fl6]). 
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First of all, we consider the situation that to is so small that (^^j = w" 1 has solu- 
tions. a n has numerous kinks (discontinuities), from blunt ones to sharp ones, according to 
Eqs. ([[5]), ([16]) and ([19]). Let us consider kinks which satisfy > i> ^ ^max)- 

From now on, we call such kinks "big" kinks. The interval between two kinks with sharpness 
0(ip) is roughly given by \^^f) ■ Thus the typical interval of big kinks is about u -1 . 
First, we divide the integration range of Eq. ([6]) into short intervals of length ~ u^ 1 around 
each big kink as 

4(*o = E J +#), ( Bl ) 
i 

where the integer / labels each big kink and P +l {k) denotes the contribution to I\{k) from 
the Z-th interval. Each interval contains one big kink and numerous "small" kinks, which 
satisfy (^^jf) < 'ip Vw*)- Let us assume that in the Z-th interval a n can be 

decomposed as 

a H (u) = af(u) + 5a?(u). (B2) 

where d' l (u) denotes the smooth function (except one big kink) which we can get after 
averaging contributions of small kinks to a H , and Saf is the contribution of small kinks. 8a'{ 
discontinuously jumps at each small kink and the width of the jump is ~ ^ 1//2 . Its average 
vanishes ((Saf) = 0) since the jump at each kink takes random values. Then we get 

a i (u) = a\(u) + 5a\(u), (B3) 

after the integration of Eq. ( 1B2I) . Then, 

I+ i(k) — / dua! l t exp(iu(u—n-a.i(u)—n-6a.i(u))/2)+ / duSa'l exp(iu(u— n-a;(it)— n-5aj(u))/2). 
Ji Ji 

(B4) 

Here the integral is performed over the Z-th interval. 

We are interested in ensemble averages of products of two of P + , for example, (l-Z^I 2 ). It 
contains mean squares of J+j's and cross terms of different J+j's. First, we evaluate mean 
squares of ilj's. (\I]_i\ 2 ) contains mean squares of the first and second terms in Eq. (IB4I) . 
and the averages of the cross terms between them. The latter vanish since (5af) = 0. In 
order to estimate the mean square of the first term in Eq. ( 1B4I) . we approximate it as 

/ duaf exp(z'u;(u-n-a,(u)-n-6a;(u))/2) ~ ^ a? (u,. a )e <Ij(Ui "'- n - a(u,i - ))/2 e- <wn - 5a(Ui "' )/2 Auj ;a . 

^ l a 

(B5) 
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Here we write the position of the a-th small kink in the l-th interval as ui- a , and Aui- a = 
ui-a+i — ui- a . We assume that each interval between two small kinks is so short that the 
integrand can be regarded as constant. We want to evaluate the mean square of this quantity, 



a,b 

(B6) 

Here we separate the average related to af , a\ and that related to Saf, 8a\, assuming that 
there is no correlation between small kinks and big kinks. In order to evaluate the second 
parenthesis in Eq. ( 1B6I) . we decompose 8a\ as 



Sat = E F i 



(k)i. 



u 



(B7) 



where F^ l {u) is the contribution of kinks of sharpness ip ~ ipk, so it has discontinuities at 



intervals ~ ( "^fc-^r (tpk) ) and between two of them its absolute value ~ ip k ' . Then 



1/2 



K|£a<(ii,. a ) - 5aK^ ;6 )| 2 )] 1/2 ~ [(\Sa}( Ul , a ) - 6a\(u l]b )f)} 



2\ll/2 



y k S / 



1/2 



1 1/2 



k s 

1/2 

1/2 



1 1/2 



'/'max / J AT 

dW- 1 ( (V) 



fdN 



\ dip 

-1/2 



\Uha ~ Ui- b 



Ul-a - Ui- b \ 

1/2 

\ui- a ~ Uhb 

1/2 



1/2 



(B8) 



To proceed from RHS of the first line to the second line, we regard F^ l (u) as constant in 
the interval between two small kinks. Here u k s denotes the position of s-th discontinuity of 
Fi(u). Fi k ^ l (u k s ) can be thought of as a probability variable whose average is and whose 

, (uK)F, ^(u k ' sl )) = unless k — k',s — s', assuming 



variance is ~ ipk- We can set (iy (u k s )F t 
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that different kinks are not correlated. This enables the second line to be simplified to the 
third line. Then we substitute ipk into (F^ 1 (uf s )) 2 , and fV'fc^rC^fc)) u i )S +i~ u i,s- The 

third factor in the forth line represents the number of small kinks in the interval (ui ;a ,ui-b). 
When we proceed from the fifth line to the sixth line, we changed the sum to the integral 
/ d(]nip) = J d^i)- 1 . Using Vw^Vw) = w and \u l a - u[\ < u' 1 , we find 

«|^K a ) - 5 ai (u, b )\ 2 ))^ < « a;" 1 . (B9) 

Therefore, wn • (5ai(ui ;a ) — 5&i(ui-b)) is much less than unity and (e lwn "^ a ^ Ui i" _,5a! '-" i;a ^' 2 ) is 
~ 1. Then Eq. (IB6I) is written as 

J2(at(ui-,a)af(u l;b )e™( Ul -- a -^ =11 dud'i(u)e iuj{u -™ {u))/2 

a,b \ Jl 

(BIO) 

The result is same as that derived without the contribution from small kinks. Then, the 
RHS of Eq. ( 1B10I) can be calculated as Eq. ( 1261) . and its magnitude is 

' ,l/2\ 2 

V) • (Bll) 

where ipi denotes the sharpness of Z-th big kink. 

In order to estimate the mean square of the second term in Eq. (IB4I) . we approximate it 

as 



duSaf exp(iw(«-n.a I («)-n. < ya,(ti))/2) ~ £ £ exp^t)« +1 -<), (B12) 

k s 

where 5f s = oo(uf s — n • ai(uf s ) — n ■ <Saj(wf s ))/2. The mean square of Eq. ( 1B12I) is 

EE F / fe)i (O ex pK s )« 

V ft s 
ft s 



+1 - U l,s) 



n,s+i u i,s) 



dN , , ,V 2 / 1 , / , dN 



£(»>) 



, -i 



ft 



^(tx) ) w" 1 . (B13) 
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Remembering ^^(^1 = k> and ip max < ^z, we find (lB13|) < flBll[) . Eventually, the 



mean square of I\ l is roughly estimated as (IBlip . In other words, in each interval around 
each big kink, it is sufficient to consider only the isolated big kink, while neglecting small 
kinks. 

Next, we consider the cross terms of different P + t s, such as (71 m ). This should vanish, 
and we can explicitly check this by straightforward calculation. To do so, we divide I % + 1 and 
71 m as Eq. (1B4|) and evaluate the mean squares and the averages of the cross terms of the 
two term, using above approximations, such as Eqs. ( 1B5I) and ( 1B12I) . 

As a result, the mean square of P + can be evaluated by summing up Eq. ( 1B11I) for each 
/. Then we obtain 

<|7;| 2 >~^^- 2 



« n 



dip^-(ip)ipuj 2 x L 
dip 



dN , -2 / 

~ fmax-^l^max) X 1p max UJ X L, (B14) 

where L denotes the integration range of P + . This implies that the greatest contribution 
to (|7^| 2 ) comes from kinks which satisfy ip ~ V'max- Such kinks dominantly contribute to 
GWs with frequency ~ uj. 



So far we have discussed the case where ipkjj^(ipk) ~ oj has a solution. However, if u is so 
large that V'fe^'(V'fe) ^ ^ is satisfied for arbitrary values of sharpness, all kinks are thought 
of as "big kinks" . Therefore, the contribution from each interval of length ~ u around each 
kink becomes ( 1B11I) . That from regions far from any kinks is exponentially small when 
uj — > oo. Eventually, 

~ ^max3T(^max) X 1p max UJ~ 2 X L, (B15) 

dip 

where ^ max denotes the value of ip at which ^r(^) has a peak. This implies that kinks which 
satisfy ip ~ f/Wx dominantly contribute to (|71| 2 )> and GWs of frequency ~ uj. Thus we 
have proved the validity of Eq. 
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Appendix C 



Here we discuss a subtlety related to validity to use the distribution function of kinks 
[Eqs. (I15p . ( fT5"j) and (|19p ]. These formulae are derived without considering gravitational 
backreaction. The distribution may be altered if such an effect is taken into account. It may 
be necessary to define the residual lifetime for blunt kinks and set lower cutoff of sharpness. 
It is difficult to clarify how we should take into account this effect at this moment. However, 
at least we can find a crude condition which must be satisfied regardless of the detail of 
backreaction; the energy of GWs emitted from strings must be less than the string energy. 
This condition is expressed as 

,dP 



dJ , 
duj' 



xt< yt. (CI) 



tot 



LHS represents the energy emitted from kinks on one infinite string in a Hubble horizon 
per Hubble time, and RHS denotes the energy of one infinite string in a Hubble horizon. 
(Note that this is only a necessary condition that the backreaction does not affect kink 
distribution.) 

First, let us assume that strings emerged in the radiation era. In the radiation era, the 
condition ( 1C1I) is satisfied for 

t < (10G/i) 3/£; % ~ (10Gfi)~ 7 -% ( <=> T > (10G/i)- 3/2 ^T, ~ (10G/i) 3 - 9 T,). (C2) 

(E = SC~ P, E r ^ -0.38, E m ~ -0.4.) For Eq. (j£5J) to be satisfied in the whole radiation 

era, 

T* < (10G^ 2E "T eg ~ (10G/i)- 3 - 9 T eg . (C3) 

If we take Gfi ~ 10~ 7 , this becomes T* < 10 14 GeV. In the matter era, the condition (1 1 j) 
is written as 

(rp \ —2E r /Em / rp \ 1.9 

jr) *eq~ (lOG^) 75 t cq . (C4) 

The condition that the backreaction is not problematic in the matter era also leads to (IC3I) . 
Eventually, if we assume that kinks had not appeared until friction domination ended or 
strings emerged at low temperature at which friction can be neglected, the gravitational 
backreaction is not important. 

Next, let us assume that strings were born at the end of inflation. It is easy to see 



that ()Cip is satisfied in the first matter era, using Eq. (145 p . After the first matter era, the 
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situation depends on whether the reheating temperature exceeds T c [Eq. (|T4l) ] or not. If 
T r < T c , there is no period when the friction works and kinks produced in the first matter 
era survive. The peak of ^xjltot consists of contribution from kinks produced around the 
turning point from the first matter era to the radiation era. Therefore, the above discussion 
applies and the condition (1C1I) is satisfied all the time. On the other hand, in the case of 
T r > T c , the friction becomes problematic in the early stage of the radiation dominated 
era. If all kinks disappear in this stage and kinks restart to emerge at the end of friction- 
domination, the condition (1C1 [) is never violated as discussed above. In the opposite case, 
where all kinks survive the friction-dominated era, (101 [) is not guaranteed. In such a case, 
the largest contribution to J dco^\ to t comes from kinks produced around the end of the first 
matter era. The condition ( 1C1I) is satisfied if 



For Gfi = 10 7 , this leads to T r /M p \ < 10 4 . Therefore the gravitational backreaction might 
be able to be neglected unless the reheating temperature is so high. 
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